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We study the effect of a single impurity on the transport properties of a one dimensional quantum 
liquid highly excited away from its ground state by a sudden quench of the bulk interaction. In 
particular we compute the time dependent dc current to leading order in the impurity potential, 
using bosonization, and starting both from the limit of a uniform system, and from the limit of 
two decoupled semi-infinite systems. Our results reveal that the nonequilibrium excitation of bulk 
modes induced by the global quench has important effects on the conductor-insulator quantum 
phase transition, turning it into a crossover for any small quench amplitude, and destroying the 
exact duality between the conducting and insulating fixed points. In addition, the current displays 
a faster decay towards the steady-state as compared to the equilibrium case, a signature of quench- 
induced decoherence. 


I. INTRODUCTION 

Transport phenomena in strongly correlated quantum 
systems are typically sensitive to interactions, inhomo¬ 
geneities and low-dimensionality and often reveal intrigu¬ 
ing and unexpected effects. A prominent example is the 
impact of a single potential barrier in a one dimensional 
interacting electron liquid. While the clean system would 
display, despite interactions, ideal conductance quantiza¬ 
tion in the presence of metallic leadsA—, the impurity at 
zero temperature is capable of turning the perfect con¬ 
ductor into an insulating link, giving rise to power law 
corrections to transport coefficients at finite temperature 
or finite voltage bias^—. The physics of a single impurity 
in a Luttinger liquid has attracted enormous attention 
in the past decadesA2£ and still finds beautiful applica¬ 
tions in systems such as the edges of a quantum spin 
Hall insulato r 11 : 12 or in electronic quantum circuits^ - — . 
Traditional condensed matter transport settings involve 
a system initially in thermal equilibrium which is per¬ 
turbed away from it by the action of some external field 
conjugate to a conserved current. Within linear response 
theory one then obtains transport coefficients, which give 
information on the structure of the underlying ground- 
state and its low-energy excitations. 

Recently, experimental advances in controlling and 
probing ultra-cold atomic gases has offered a new plat¬ 
form to study nonequilibrium time-dependent phenom¬ 
ena in a fully tunable settin g 16 : 17 . First generation of 
experiments studied the dynamics induced by rapidly 
changing in time some system parameter—in an other¬ 
wise isolated system, a so called quantum quench. More 
recently, the experimental focus shifted towards realizing 
genuine transport experiments with cold atomsSL- 22 . It is 
important to realize that these systems are very well iso¬ 
lated from their environment and hence intrinsically out 
of equilibrium, therefore the standard condensed matter 
idealizations do not necessarily apply. 

Motivated by these experimental developments, in this 


paper we investigate the transport properties of nonequi¬ 
librium quantum many body states which are thermally 
isolated and highly excited above their ground state. 
While for generic, non-integrable and ergodic, quantum 
many-body systems one may expect the excitation en¬ 
ergy to be effectively converted into temperature at long 
time s 23 : 24 , and hence to induce decoherence, the situation 
is less clear at intermediate time scales where long-lived 
prethermal states may emerge^—, especially paramet¬ 
rically close to an integrable point. The idea we pursue 
in this work is to use transport as a probe to unveil the 
structure of these prethermal excited states, to under¬ 
stand the relevant excitations and whether non-trivial 
quantum phenomena survive at these high-energies, and 
if so, in which form. It should be noted that there has 
been a recent spur of interest towards understanding sim¬ 
ilar dynamical quantum correlations in isolated many 
body systems excited after quantum quenches^ - — . We 
have considered an example along these lines in a re¬ 
cent work on a time-dependent orthogonality catastro¬ 
phe problem^ where we introduced a novel dynamical 
Loschmidt echo, encoding the response of a highly ex¬ 
cited state to a local perturbation. Here we study a sim¬ 
ilar problem from the point of view of transport. 

Specifically, in this paper we compute the time- 
dependent current of a one dimensional system of spinless 
electrons described by the Luttinger model, which is ex¬ 
cited by a sudden change of the two-particle interaction 
and a simultaneous switching-on of a local scattering po¬ 
tential. These two perturbations have very different ef¬ 
fects, the former injecting extensive energy into the bulk 
modes, the latter creating a non-linear channel for local 
dissipation. We discuss the physics both in the limit of a 
uniform system, where the local potential induces a weak 
backscattering term, and in the limit of two decoupled 
semi-infinite systems where we switch on a local tunnel¬ 
ing. In both cases, using bosonization, we formulate the 
problem in terms of a boundary sine-Gordon model in 
a nonequilibrium transient environment. Using pertur¬ 
bative approaches we discuss the fate of the conducting- 
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insulating zero temperature quantum phase transition in 
the presence of a nonequilibrium bulk excitation, thus 
complementing and expanding previous results 37,38 . 

The results for the steady state current reveal the 
emergence of a novel energy scale associated with a 
quench-induced decoherence effect, which turns the sharp 
equilibrium transition into a smooth crossover. Interest¬ 
ingly, the very same energy scale has been shown to play a 
key role in the transient orthogonality catastrophe prob¬ 
lem, cutting off the renormalization group flow of the 
backscattering and resulting in an exponential decay of 
the Loschmidt echo^. 

Our analysis, although perturbative in nature, ulti¬ 
mately suggests that the steady state impurity problem 
in a quenched Luttinger model has very different proper¬ 
ties than in equilibrium at zero temperature. In particu¬ 
lar the quench acts as a relevant perturbation on both 
sides of the equilibrium transition, even when the lo¬ 
cal non-linearity would be irrelevant at the tree level, 
generating an effective decoherence which drives the sys¬ 
tem away from the uniform and open-chain fixed points. 
While this behavior is reminiscent of an effective tem¬ 
perature, we will see that this analogy is only qualitative 
rather than quantitative, as the scaling of physical quan¬ 
tities in the steady state with respect to this emergent 
energy scale do not show signature of non-trivial power- 
laws like the ones encountered at non-zero temperatures. 
At the same time we find that deviations from the asymp¬ 
totic steady state regime, due to transient effects at finite 
time, display nonequilibrium power laws with character¬ 
istic Luttinger liquid exponents which may be a signature 
of a sort of Luttinger liquid universality out of equilib¬ 
rium, as recently discussed 29 ! 38 ' 39 . 

Note that the problem of an impurity in a quenched 
Lut ting er model has been recently addressed in Refs. [37l 
andl38l. In the former, the dynamics of the current in the 
tunneling regime was investigated after a global quench of 
the bulk interaction and a faster decay to the stationary 
state was found. In the latter, results for the tempera¬ 
ture scaling of the conductance in the steady state after 
the quench have been obtained using different approaches 
such as bosonization and functional RG. Although we do 
not address temperature dependence of transport coef¬ 
ficients in this paper, and we work always at zero tem¬ 
perature, it is useful to comment on the relation between 
the results of Ref. [H and the physical picture emerging 
from this work. We will do this in the discussion section. 

Finally, it is worth mentioning that the question of 
nonequilibrium effects on local quantum criticality has 
been addressed recently also in the context of driven 
quantum systems^ and specifically for single impurity in 
Luttinger liquids in the context of noisy driven Joseph- 
son junction circuits'- 1 -' or tunneling in biased quantum 
wires - 2 using the out-of-equilibrium bosonization frarne- 
wor k 43,44 . The picture emerging in these cases, where 
the nonequilibrium perturbation induces a decoherence 
mechanism which eventually cuts the quantum critical 
power-law behavior, is consistent with our results for the 


isolated quenched problem. 

The paper is organized as follows. In section |IT] we in¬ 
troduce the system, the nonequilibrium protocol, and a 
derivation of an expression for the time dependent cur¬ 
rent in terms of a bosonic Green’s function. We then eval¬ 
uate the current in the weak-backscattering limit in sec- 
tion lll Al Such a perturbative approach eventually breaks 
down in a certain region of parameter space and this will 
lead us in section m to formulate the problem in the 
limit of two disconnected ID systems coupled by a local 
tunneling term. We will compute perturbatively the tun¬ 
neling current in section llll Al and thus obtain a complete 
picture of the problem. 

II. TRANSPORT IN THE 
WEAK-BACKSCATTERING REGIME 

We start by discussing the model and the nonequilib¬ 
rium set up. We consider a one dimensional system of 
interacting spinless fermions described by the Tomonaga- 
Luttinger (TL) model. We assume the system to be ini¬ 
tially (at time t < 0) in the ground-state |\f r o) of the TL 
Hamiltonian 

H 0 = H {iee + y J dx (p 2 L (x) + p R (x)) 

+3° J dxp L (x)p R (x) (1) 

where Uf ree = v F J2 a =L,R I dx^ a id x i’ a with s L = 
—sii = 1 , p a = il>aipa is the fermion density with given 
chirality a = L/R , and 34 ( 3 !)) are the strengths of the 
intra (inter)-branch scattering processes. 

At time t = 0 we assume a sudden change (quench) in 
the value of the bulk interactions to 32 , 34 , so that the 
system will evolve with the Hamiltonian 

H = L/free + y J dx (Pl(x) + p R {x)) 

+32 J dxp L (x)p R {x ) (2) 

This global quench injects extensive energy into the sys¬ 
tem and triggers a nonequilibrium occupation of the bulk 
modes. Dynamics arising simply due to this quench has 
been discussed extensively in the literature^—. In ad¬ 
dition to the bulk excitation, we assume that a static 
impurity potential Vj oc has also been suddenly switched 
on at time t = 0. This static impurity can induce an 
intra-branch (I?, L ) — > ( R , L) as well as inter-branch 
( R,L ) -* (L,R) scattering of fermions, that we write 
as the sum of two contributions 

Moc = Vb (4 (0)^L (o) + L o r) + (3) 

+^bs (4(0)V’h(0) + h.c.'j . (4) 

the first term representing the forward scattering while 
the latter the backward scattering. As a result the wave 
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function at time t is |\t(i)) = exp (—iH+t) l’F(O)}, with 
H + = H + V loc . 

We employ bosonization to describe the system, thus 
introducing the bosonic fields 4>(x) and 6(x) describing 
collective density and current excitations respectively, 

^ L/R ( x )r,^ = e i ^ ± ^ (5) 

\J I'KOi 

Pl/r{x) = -d x (<f>(x) ± 6(x)) /2tt (6 ) 

The Hamiltonians before and after the quench in terms 
of these bosonic fields are 


Hn = 


Zio 

2tt 


dx 


K 0 (d x 9(x)) + — (d x <t>(x)) 
Ao 


( 7 ) 


and 


H + = 


27T 


dx 


I< ( d x 6(x)) + — {d x <j>{x)) 

K 


+ ^loc (8) 


Within bosonization, the impurity potential in Eq. © 
can be written as 


Vioc = g is d x <))(0 ) + g hs cos 2</>(0) (9) 


where the effective impurity couplings read respectively^ 


9 fs — 


Vh 

7r 



7ra 


( 10 ) 


The Luttinger parameter I\ o, K and the sound velocities 
Uq , u are related to the Fermi velocity and the interaction 
parameters as 


As usual, the current has two contributions, J(x,t) = 
Jd{x,t) + J p (x,t), the diamagnetic ( Jd ) and the para¬ 
magnetic (J p ) one. The former is given by (hereafter 
h= l) 


g 2 uK 

Jd{x,t) = - A(x,t) =—DA(x,t) (16) 

7T 

where we have introduced the diamagnetic term D = 
e 2 uK/TT. The latter is 

Pll r\ P 

J p (x,t) = - d x 0(x,t) = -d t tp(x,t) (17) 

7T 7T 

Finally, we notice that the total current can be equiva¬ 
lently obtained from the continuity equation 

d t p(x, t) -(- d x J(x,t) = 0 

by computing the time derivative of the density p(x, t) = 
-(e/7r )d x (f>(x,t) 

1 P 

d t p(x,t) = - [H(A),d x cj>(x,t)] 

7T 

using the LL Hamiltonian with the minimal substitution. 
The result of this calculation obviously recovers the ex¬ 
pression for the current J(x,t) = Jd(x,t) + J p (x,t). 

We can now compute the average current within linear 
response theory, by doing perturbation theory in A(x , t). 
Since the diamagnetic part of the current Jd(x,t) is al¬ 
ready linear in the vector potential, we only need to take 
into account the paramagnetic contribution which gives 


u 0 = v F \J (l + gl/2itv F f - (gl/2-Kv F ) 2 (11) 

K _ l l + gj/2Tw F -g°/2nv F 
0 V l + g°J‘2nv F +g° 2 /2irv F 1 j 

with similar relations holding for u, K as a function of 
< 72 , < 74 - In order to preserve Galilean invariance we choose 
the values of < 72,54 such that uqKq = uK. This amounts 
to simply requiring that <74 — 52 = 54 — 5° ■ 

In order to probe transport through the system we 
will study linear response to a weak electric field E(x, t) 
applied after the quench, 

E{x,t) — —d t A(x,t) (13) 

with A(x, t) being the vector potential. In the presence 
of the electric field, the LL Hamiltonian is modified ac¬ 
cording to the minimal substitution which amounts to 
the shift 


{Jp{ x, t)} 



dx'x$j{xt;x't')A(x',t') (18) 


where we have defined the retarded current-current cor¬ 
relation function 

XjjOrf; x't') = —iO(t - t'){ [ J p {x , t),J p (x', t')}) 

= -i-—r0(t - t'){ [dt4>{x,t),d t '(l){x\t')]) (19) 

7T 

and extended the time integral up to infinity. We stress 
here that the average in Eq. © is taken over the initial 
density matrix (ground state of H 0 ) while the operators 
are evolved with the Hamiltonian H + . The resulting cor¬ 
relator is not the usual equilibrium one and in particular 
it is not time-translational invariant. This reflects the ef¬ 
fect of the global and local quantum quenches that have 
been performed on the system at time t = 0. Finally we 
have for the current 


d x 0{x) —> d x 6{x) — eA(x, t) (14) 


(J(x,t)) = —DA(x,t) + {J p (x,t)) (20) 


The current operator can be obtained as the functional 
derivative of H{A ) with respect to the vector potential 


SH 


(15) 


We now use an important identity relating the current- 
current correlation © to the retarded Green’s function 
of the field <j>(x, t) defined as 

Q R (xt\ x't’) = —i0(t - t')( [4>{x, t), 4>(x',t')]) (21) 


J(x, t) 


5A(x, t) 














4 


We have 

P 2 

Xjj(xt', x't') = -^d t 'd t Q R {xt\ x't') — 2 DS(x — x')8(t — t') 

7 T Z 

( 22 ) 

Substituting this result into Eq. m and using 
Eq. (HU), we find that the last term cancels the diamag¬ 
netic term exactly, and after an integration by parts we 
obtain 

(J(x, t)) = — -^2 J dt' j dx' d t Q R {xt] x't') E{x't') (23) 

This is in principle an exact result for the time depen¬ 
dent current through a finite-size LL and we can now 
further specify the electric field profile. We assume a 
sudden switching of the field, E(x,t) = 9(t)E(x) and in 
addition, since we are ultimately interested in the cur¬ 
rent due to the impurity ©, we assume that the po¬ 
tential drop occurs around x = 0, V{x) = V6(— x), 
so that the electric field is effectively a delta-function 
E(x,t > 0) = — d x V(x) = V5{x). Then, the current at 
x = 0 , (J( 0 ,f)} = I(t) can be written only in terms of 
the exact retarded Green’s function of the local field at 
the impurity site, 

G R {t,t') = -iO{t-t')m)A{t')}) (24) 

with cj)(t) = </>( 0 ,f), and it reads 

I(t) = f dt'd t G R (t,t') (25) 

ft Jo 

This equation, which gives the time-dependent current in 
terms of an exact dynamical correlator of the local field, is 
the main result of this section. It generalizes to the time- 
dependent quench problem the well known equilibrium 
result^. From this we can compute the current for a 
pure LL, using the expression for the bare retarded local 


Green’s function 

G R (t > t') = — K arctan A ( t — t') (26) 

where A = u/a is an ultra-violet cut-off. We obtain 
e 2 KV 

Io(t) = - 5 —arctan At (27) 

7 T Z 

which approaches in the long time limit Io ss /V = 
e 2 K!2Tt = e 2 K/h as expected. In the next section we 
will evaluate the corrections to the time dependent cur¬ 
rent due to the local potential. The forward scattering 
term in Eq. ([5]) turns out to not contribute to the current, 
as we discuss explicitly in appendix E therefore in the 
next section we will focus our attention on the backward 
scattering term which is responsible for the interesting 
physical effects with V\ oc = gbs cos( 2 ^( 0 )). 

A. Weak-Backscattering Corrections to Linear 
Conductance 

We are now in the position to evaluate the weak back- 
scattering correction to the conductance. We just need 
to evaluate the local Green’s function for the field <f>(t ) 

G a b(t,t') = (28) 

where a, b = ± are Keldysh indices. We have defined the 
various Green’s functions and the relations between them 
in appendix [Aj To lowest order in the backscattering gbs 
we obtain^ 

8G a b(t,t') = f dtidt 2 y^ j (cd)F cd (t 1 ,t2)Kb( t i,h) 

cd 

(29) 

where we define 5G a b = G ab - G 0 ,ab, 

F c d(ti,t 2 ) = —i(cos 2(j) c {t\) cos 2(f) d (t2)) 

and 


Kb(G,t 2 ) 


— Go,ac(t, l)Go,cb(l, t') + Go,ad(t, 2)Go,db(2, t') — Go,ac(t, l)Go,db(2, t') — Go, a d{t, 2)Go,cb(l, t') (30) 


The retarded component is given by 


G R {t,t') = ~Y.bG ab {t,t'), 

ab 


lit) - Ioit) 


(31) 


8I(t) = —C f dfi<9 t G^(f,fi) f dt 2 F R it 1 ,t 2 )$iti,t 2 ) 
Jo Jo 

(32) 

with G = FJlV- g ? s an d 



Substituting the resulting expression into Eq. (l25l) . and 
after some algebra we obtain for the current Slit) = 


<f>ih,t 2 )=K 


dr arctan At = KAt\ arctan Ati + 
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FIG. 1. Parameter space for the weak-backscattering prob¬ 
lem. In the white area (region I) corresponding to A' neq > 1, 
perturbation theory is well behaved, while the shaded area 
(region II) corresponds to the regime where backscattering is 
a relevant perturbation in the steady state, and the pertur¬ 
bative expansion breaks down. The dashed line is a reference 
only to the validity of the strong coupling (weak-tunneling) 
expansion (see further below). 
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FIG. 2. Backscattering corrections to steady state current 
(top panel) as a function of Ko and for different values of 
K. A finite quench amplitude restores transport with devia¬ 
tions from Ko = K (arrows) vanishing quadratically for small 
quenches, 8I SS /V ~ —g 2 B (K — Ko) 2 (see bottom panel, log-log 
scale). Parameters: A = 1, </ba = 0.1 


—KAt 2 arctan At 2 —— log 

The retarded correlator entering the expression for the 
current is found to be 


fl + AHj\ 
\ 1 + A 2 t| / 


F«( h > t 2 ) = _ ^(2A-arctanA(f 1 - f 2 )) /ft ^ 2) 


(l + A2 (h-hf) 

where we have introduced 


(33) 


f(h,t 2 ) 


(l + 4A 2 t\) (l +4A 2 tf) 
(1 + A 2 (t\ + t 2 ) 2 ) 


(34) 


and the exponents 



We can further simplify Eq. (1321) by noticing that for 

Mi > 1/ A 


d t Go(t - h) ~ - 7 xK8(t - ti) 

Then we can write the correction to the time dependent 
current as 


SI(t) =~C dr 


sin [2 K arctan Ar] 
(1 + A 2 r 2 ) Kneq 




T) $(t,t-T) 

(36) 


with C = nKC = 4e 2 A'Kg 2 s /7r. In the next section we 
will discuss the behavior of this quantity as a function of 
time and quench amplitude. 


B. Discussion 


We start by discussing the equilibrium zero tempera¬ 
ture case corresponding to K 0 = K. In this case we have 
K tI = 0 and A' neq = K and the integral in Eq. (1551) sim¬ 
plifies. We obtain for the transient current at long times 
(At > 1) 

(37) 

i.e., for K > 1 the weak-backscattering correction to the 
current vanishes in the long time limit as a power law, the 
junction is perfectly conducting and the backscattering 
is irrelevant, this is the well known result from Kane and 
Fisher—. As soon as the quench amplitude is non-zero, 
K Ko, we find a number of interesting differences. 
We first focus on the long time steady state value of the 
backscattering current, SI ss = I ss — Iq ss , this is 


SI ss /V = -^K 2 gl s G 0 / dr 


sin (2 AT arctan Ar) 

(1 +A 2 r 2 )^ neq 


(38) 

Go = e 2 /27 rh = e 2 //i the quantum of conductance, which 
we set to one from now on, Go = 1. We first notice 
that the integral is well defined as long as AT neq > 1, as 
expected from perturbative RG calculations^, since only 
in this regime the local backscattering is an irrelevant 
perturbation for the steady state. We plot in figure |T| 
the region of parameters where this condition is satisfied, 
and our calculations are therefore controlled. In figure [5] 
we plot the backscattering correction to the steady state 
current for different values of K in region I, as a function 
of Kq. Quite interestingly, the correction to the steady 
state current is different from zero at any finite quench 
amplitude AT ^ Kq , in other words the systems deviates 






















from the perfect conduction limit. For a small quench 
amplitude one finds that the steady state correction to 
the conductance is proportional to ( I\ — Kq) 2 , 

SI ss /V ~ -g 2 s ( K - Kq ) 2 (39) 

We now discuss the transient behavior and the ap¬ 
proach to the steady state. To this extent we evaluate 
numerically the integral in Eq. (1361) and plot the result in 
figure [3] We find that the current decays to the steady 
state in a power law fashion, 

61{t) -5 I bb ~ i (40) 

with an exponent a(K, Kq) whose dependence on the 
quench parameters is shown in the bottom panel of fig¬ 
ure [3] as a function of Kq at fixed K = 1.25 in the 
region I. We notice that the exponent a behaves non- 
monotonically with Kq and reaches a minimum value for 
K = Kq , i.e., the decay of the current is stronger in 
the presence of a finite quench amplitude. We can get 
an analytical understanding of this by looking at the in¬ 
tegral expression for 51 it), Eqn. (13(1 . Using the fact 
that for large £ we have f{t,t — r) ~ 1 + 0(l/t) while 
$(£,£ — r) ~ t + 0(l/t), we obtain a = 2 (K neq — 1). 
We notice that this is indeed confirmed by the numeri¬ 
cal evaluation of the integral, although for larger values 
of the bulk quench, deviations start to appear that we 
attribute to the finite time resolution of the numerical 
integration. 

We have analyzed only a few values of I\ while tuning 
K 0 , yet we expect this behavior to hold throughout the 
region I in figure [T] We however expect the perturbative 
expansion to eventually break down as we approach the 
regime K neq = 1 and enter region II, where the backscat- 
tering becomes a relevant perturbation, i.e., its strength 
grows under RG. In equilibrium, a complementary ap¬ 
proach to explore this strong coupling regime amounts 
to starting from the limit of two decoupled semi-infinite 
systems and switching on a local tunneling. In the next 
section we will suitably generalize this approach to the 
quench case and discuss to what extent it can be used 
to extract useful information about the strong coupling 
regime. 


III. TRANSPORT IN THE TUNNELING 
REGIME 


We now consider a different quench protocol which 
will allow us to access the strong coupling regime where 
the back-scattering is relevant. Let us suppose at £ = 
0 we have two disconnected and identical semi-infinite 
TL models describing ID interacting spinless fermions. 
By definition, each lead contains left and right moving 
fermions interacting in the bulk through intra-branch 
(<?2) and inter-branch (g$) scattering processes, as in the 
previous case. A major difference however arises due to 



_k 0 = o. 25| : ny i 1 i 1 i 1 i 1 i 1 i r 



FIG. 3. Power law decay in time of the weak backscatter- 
ing correction to the current, for fixed K = 1.25 and differ¬ 
ent values of I\o■ We notice the exponent a depends non- 
monotonically on Kq and takes a minimum value in the equi¬ 
librium case K = Kq (see bottom panel). 


the presence of a sharp edge in each lead, which imposes 
open boundary conditions for the fermionic field. As a re¬ 
sult, and differently from the translational invariant case 
of the previous section, the two fermionic species are not 
independent of each other. This suggests an equivalent 
and more convenient representation of each lead * = 1,2 
in terms of a single chiral fermionic field ipi(x) defined 
on an infinite system^ (i.e., obeying periodic boundary 
conditions). In terms of this field, the initial Hamiltonian 
in each lead becomes 


Hq i — L/fr, 


/ dxpi(x)pi(x) + ^ dxpi(x) Pi(-x) 


(41) 

where Ffree =vfJ dxipjidxipi, while (x) = i/>J(x)i/>i (x), 
and one should notice the non-local g% interaction that 
results from the single chiral field representation. 

We prepare the system in the ground state | To) of 
H 0 = H 0i and then, for £ > 0, we evolve the system 
with a different Hamiltonian H 


H= Y, H i+Y, Vi N i + Hr (42) 

i=l,2 i=l,2 


where we have (a) quenched the bulk interactions, 


Hi = H fr , 


54 

2 


dx pi{x)pi{x) 


92 

2 


dxpi(x) pi ( x) 


( 43 ) 

(b) switched on a tunneling term coupling the two wires 
at the edge 


H t =Xq (v4(0)iM0) + Ac.) (44) 

and finally (c) switched on a bias voltage jii .2 = ±eU/2 
that couples to the charge imbalance. We are interested 
in computing the time dependent current 

/(£) = (£)|/|tf (£)> = (I'd e iHt Ie ~ lHt |* 0 > (45) 
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where the current operator is 


as well as the electron density in the lead i as 


I = ie[H, (TVi — jV 2 )] = — 2ieA 0 ( \p\(0)ip2{0 ) — h.c^j 

_( 46 ) 

In order to proceed, it is useful to perform a unitary 
transformation to eliminate the bias voltage in H , which 
amounts to going to the rotating frame defined by an 
operator fil(i) = exp (i JT fa iVj t). Inserting II into the 
expression for the current (1451) we obtain 

m = m)\i(t)\m) m 

where the new state 

| fit (*)!*(*)) (48) 

satisfies the Schrodinger equation 

id t \*(t))=H(t)\*(t)) (49) 

with a rotated-frame Hamiltonian 

H(t) = (idttf) n + = J2 H i + H T {t) (50) 

i 

While the bias voltage has disappeared, the tunneling 
has become explicitly time-dependent 

H T (t) = tf(t)H T Sl(t) = A 0 (e iey Vi t (0)V’2(0) + /i.c.) 

(51) 

where we have used the fact that 

& (t)i/ji(x)tt(t) = (52) 

The current operator in (1451) also acquires an explicit 
time dependence 


Pi(x) = d x ipi(x)/2n (57) 


The Hamiltonian before and after the global quench of 
the interactions (Eqns [44] and [43jl become now quadratic 
in this bosonic field and can be easily diagonalized with 
a Bogoliubov rotation (see appendix ICl). For example we 
have for 

H °' = f (~k + Ka )S t (d - v,ix))2 + 

/ ^ d ^i( x ) 9 x<Pi(~ x ) (58) 

The resulting harmonic theory is again described in terms 
of two parameters, the sound velocity uo, u and the Lut- 
tinger parameter /to, K, which by construction are iden¬ 
tical in the two leads, and whose expressions in terms of 
the coupling constant g ®, g 4 and <72 , 34 are the same as in 
the translational invariant case (see equations (JTTJ) ). We 
rewrite them here for convenience 


uo = v F \Z(! + g\l2-KVpf - (g%/2-jrv F ) 2 


K a = 


1 + gl/2nv F - g%/2nv F 
1 + gl/2nv_ F + g%/2nv F 


(59) 

(60) 


with similar relations holding for u , K as functions of 
52 , 34 - In addition we need to write the tunneling and 
current operator in the bosonic language. To this extent 
it is useful to introduce the combinations 


^(x) 


<Pi{x) ± <£ 2 ( 2 ;) 

V2 


(61) 


I(t) = fit(t)7fi(t) = —2*eA 0 (e ie ' / Vi t (0)V’2(0) - /i.c.) 

(53) 

The time dependent state |^(t)) can be written in the 
interaction picture Oj = with respect to 

the unperturbed Hamiltonian of the uncoupled wires as 

|'I'(t)} = e - i 'Z* Hit Te- i ti dt ' HT x(t’) |^ 0 ) (54) 


We can now evaluate the time dependent current to the 
lowest order contribution in the tunneling, which gives 


= dir <*o| //(*),-Hr, j(*i) |*o) (55) 


where all the operators are explicitly time dependent as 
in previous equations and evolved with the Hamiltonian 
of the uncoupled wires, after the global quench. 

It is now useful to use bosonization, in its open¬ 
boundary formulation^ to proceed with the evaluation 
of the tunneling current. We first write the fermionic 
operator in each lead in terms of a single chiral bosonic 
mode 


^ i (*) = 4=e i ' Pi{x) 

\Ja. 


and to notice that only the combination enters both 
the tunneling and the current. Indeed we have disre¬ 
garded the Klein factors which can be shown to be unim¬ 
portant. Thus, 

Hr(t) = 2Acos (^V2ip~(t) — eVtj (62) 

I(t) = -4eAsin (y/2tp-(t) - eVtj (63) 

where we have defined <p~{t) = <p~(0,t) and A = Ao/a. 


A. Tunneling Correction To Current 


In order to proceed further, we plug Eqns. (1621631) in 
the expression for the tunneling current to obtain, after 
some simple algebra, 


/(f) = — 8eA 2 / dr sin (eW) F H (t, t — r) (64) 


where F R (t,t') is the retarded component of the correla¬ 
tor 

Fab(t,t') = -i( cos (y2tp-(t)j cos (t/2^(t'))> = 


(56) 
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FIG. 4. Parameter space for the tunneling problem. In the 
white area (region I) corresponding to Ffneq* > 1, perturba¬ 
tion theory is well behaved, while the shaded area (region 
II) corresponds to the regime where the tunneling is a rel¬ 
evant perturbation in the steady state and the perturbative 
expansion breaks down. The dashed line is a reference only to 
the validity of the weak backscattering expansion previously 
discussed. 


= _1( f.i'SZVa M g-iv'ay’s (t')j (55) 

It is convenient to express this correlator in terms of the 
fields <pi(t) in the two semi-infinite leads. We have 

F ab (t, t') = e ivia{t) e - ^ 16 ^) ( e iip2b{t,) e _ivio(t) ) 

where we have used the fact that in the absence of tun¬ 
neling the two leads are decoupled. To compute this cor¬ 
relator we need therefore the local Green’s function of the 
field tfi after a quench of the bulk interaction parameter 
in a semi-infinite chain, which we give in appendix[U] The 
result for the retarded component, F R = — Fab/% 

is found to be 




where we have introduced 




(l + 4A H\) (l + 4A 2 f|) 
(l + A 2 (f 1+ f 2 ) 2 ) 2 


Tv-dual 

^tr 


/2 


(67) 


with the exponents 


77^dual 

-^■neq 


l 

2Aq 


1 + 


K 2 

A' 2 


7>^dual 

iV tr 


l 

2Aq 



( 68 ) 

We stress that Eq. (R>71) is perturbative in the tunnel¬ 
ing but contains the bias voltage eV to all orders. In 
fact in the absence of the bulk quench, K = A'o, we 
can recover from this the well known result^ for the non¬ 
linear I-V characteristic of the wire, I(V) ~ V 2 ! K ~ l . 


In the following we will instead consider mostly the low 
bias regime, for a finite quench amplitude and to this ex¬ 
tent we will evaluate the current in the low bias regime 
and equivalently the non-linear differential conductance 
at zero bias. Note that the latter I/V |y=o = 0 in the 
absence of quench. 

After some simple manipulations, we can write the 
time dependent current at small bias voltage as 


„ „ f* rsin arctanAr) 

I(t) = 8eW / dr Kk 

Jo (1+A 2 r 2 ) 


T ) 

(69) 


We notice the result we have obtained for the current 
in the tunneling regime is dual to the one in the weak- 
backscattering case in the sense that the two currents are 
related by a transformation of the Luttinger parameters, 
before and after the quench, into 


K , K 0 


1 1 
K'~K 0 


(70) 


This is the natural generalization of the equilibrium du¬ 
ality^. Yet, as one can immediately see from the expres¬ 
sions m for Ar neq ,A'tr and (11)51) for K^ 1 , Af£. ual , the 
exponents controlling the decay of the correlation in the 
tunneling and weak-backscattering regime are not related 
by this simple duality, i.e., 

I<X y 2 1/K neq (71) 

Atr 1 £ 1/A tr (72) 

unless of course K = I\q. This will have important con¬ 
sequences as we are going to discuss in the next section. 
A simple way to understand the origin of this result is 
to notice that while for the Hamiltonian after the quench 
the duality still holds, the initial condition of the problem 
(the ground state of the Luttinger model with interaction 
parameter Ao) is not dual when written in the basis of 
the natural eigenmodes of the system after the quench 
(Luttinger model with interaction parameter K). The 
fact that the Luttinger model is integrable and therefore 
never fully loses memory of this initial condition implies 
the breakdown of duality in the stationary state. 


B. Discussion 

As in the weak-backscattering case, we start discussing 
the equilibrium zero temperature result corresponding to 
A'o = K. In this case we have K tr = 0 and A neq = A, 
and the integral in Eq. (1691) simplifies. For the transient 
current at long times (At 1) we obtain, 

r m/v ~ (73) 

i.e. for K < 1 the correction to the current vanishes 
in the long time limit as a power law, the junction is 
perfectly insulating and the tunneling is irrelevant^. We 


















FIG. 5. Steady state tunneling current as a function of Kq 
and for different values of I\ (top panel). We see that a non¬ 
zero quench amplitude results in a finite tunneling current, 
vanishing quadratically as A'o I\ (arrows), I s $/V ~ A 2 (A' — 
Kq ) 2 (see bottom panel, log-log scale). Parameters: A = 
1, A = 0.1 


FIG. 6. Power law decay in time of the weak tunneling cor¬ 
rection to the current, for fixed K = 0.9 and different values 
of Kq. We notice the exponent a depends non-monotonically 
on A'o and takes a minimum value in the equilibrium case 
K = A'o (see bottom panel). 


now consider the finite quench case, K ^ Kq, and first 
focus on the long time steady state value of the tunneling 
current which reads 


I 8S /V = 8e 2 A 2 



dr 


t sin ( j7 arctan At) 
(1 + A 2 t 2 ) k ™* 


(74) 


We notice that the integral is well defined as long as 
-^neq 1 > 1, a condition which is satisfied in the region 
of parameters plotted in figure 0 ] (region I). Interestingly, 
this region is not simply the dual of the region I for the 
weak-backscattering case, (light shaded line inside region 
I in figure 01) as it would be in the equilibrium case. This 
is a consequence of Eqns (1711721) . Therefore we conclude 
that the quenched impurity model does not display a 
full duality, contrarily to other nonequilibrium realiza¬ 
tions such as the driven noisy one^i. We will come back 
to this point in the next section. Let us now discuss 
the behavior of the steady state current. Quite inter¬ 
estingly, the steady state current is different from zero 
at any finite quench amplitude K ^ A'o, in other words 
the system deviates from the perfect insulating limit. In 
figure 0] we plot the steady state value of the tunneling 
current at fixed K and as a function of Kq. As in the 
weak-backscattering limit we find that for small quench 
amplitudes the steady state current is 

Iss/V ~ A 2 (A" — Af 0 ) 2 (75) 


We now discuss the transient behavior and the ap¬ 
proach to the steady state. We evaluate numerically the 
integral in eq. ( 1691 ) and plot the result in figure 01 We find 
that the current decays to the steady state as a power law 

m - 4 S ~ ^ (76) 

with an exponent a(K, I\q) whose dependence on the 
quench parameters is shown in figure 01 As in the weak- 


backscattering case previously discussed, we find that the 
decay of the current is in general faster for a finite quench 
amplitude, in accordance with the analysis of Ref. IffR 
In particular the exponent a reaches its minimum for 
K = Kq and behaves as a ~ 2(A'(J“q 1 — 1). 

IV. DISCUSSION 

Let us summarize the results of previous sections and 
discuss their physical consequences. We started from the 
limit of a uniform system and found that for K neq > 1, 
corresponding to region I of fig. [U perturbation theory in 
the backscattering is well behaved. While this means that 
in equilibrium, K —> Kq, the conductance approaches 
the unitary limit, a bulk quench K ^ Kq gives rise to 
a finite deviation from this perfect conducting regime. 
This is consistent with the RG analysis we developed in 
Ref. HE There we showed that although for AT neq > 1 
the backscattering is an irrelevant perturbation for the 
steady state, it nevertheless gives a sizeable effect, gener¬ 
ating an effective temperature for the local degree of free¬ 
dom, which explains, at least qualitatively, the deviation 
from the unitary limit. When K neq < 1, the backscatter¬ 
ing becomes relevant, and perturbative approaches in the 
back-scattering breaks down, and one enters the strong 
coupling regime. 

A possible approach to grasp the behavior of the sys¬ 
tem in this limit of strong backscattering is to start 
from the strong coupling fixed point, assuming the ef¬ 
fective temperature is not enough to cut the growth of 
the backscattering under RG flow, and perturb this sys¬ 
tem of two decoupled wires with the sudden switching on 
of a local tunneling. We have considered this regime in 
section Hill and found that for A'))™ 1 > 1, corresponding 
to region I in fig. 01 perturbation theory in the tunnel¬ 
ing is well behaved. While this could naively suggest the 
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FIG. 7. Steady state current versus Luttinger parameter I\q 
at fixed values of the quench amplitude K — Kq, in the tunnel¬ 
ing regime (left panel, for A = 0.1) and in the weak backscat- 
tering limit (right panel, for g^ a = 0.1). We see that the sharp 
transition between insulator and conductor driven by Kq gets 
smeared out by the global quench. 


strong coupling fixed point remains stable as in equilib¬ 
rium, we have found that a finite tunneling current ap¬ 
pears in the steady state for any finite quench amplitude 
K ^ Kq. This result is again reminiscent of a thermal 
behavior and we may speculate that an effective temper¬ 
ature would indeed be generated under RG by the local 
tunneling, in analogy with the backscattering case. 

Taken together, these results suggest that a global bulk 
quench has a dramatic effect on the problem, smear¬ 
ing out the sharp distinction between the conducting 
and the insulating phase that exists in equilibrium at 
zero temperature. To further appreciate the role played 
by the quench amplitude in the problem, it is use¬ 
ful to look again at the steady state current in both 
regimes. In figureQwe plot the weak-tunneling and weak- 
backscattering correction to the current as a function of 
K and at fixed K — Kq. At zero quench amplitude this 
would give the usual sharp transition from an ideal insu¬ 
lator for K < 1 to a perfect conductor for K > 1. How¬ 
ever as we clearly see, any small finite quench amplitude 
turns this into a smooth crossover. Finally, it is worth 
noticing that, as we mentioned earlier, the two regions 
A' neq > 1 and > 1 are not dual to each other. This 

means that there exist regions of the parameter space 
(K, Kq) where neither the backscattering nor the tunnel¬ 
ing are relevant, i.e., perturbation theory is well behaved. 
We interpret this result as a further signature of the ef¬ 
fective temperature behavior of the problem. 

These results for transport offer therefore another ex¬ 
ample of the quench-induced decoherence mechanism 
that we had previously identified in the orthogonality 
catastrophe problem^. As in that case, the relevant en¬ 
ergy scale, in the limit of small quenches, is proportional 
to the quench amplitude squared 

7 * - g 2 (K - Kq) 2 (77) 


where g in the equation above indicates any source of 
local non-linearity, either due to the backscattering po¬ 
tential in the weak impurity limit, or due to the tunnel¬ 
ing term in the limit of two decoupled chains. While the 
qualitative behavior is suggestive of an effective thermal 
behavior, it is important to stress that such an analogy 
does not fully carry over. In particular, the characteristic 
power law structure exhibited by the finite temperature 
current in the equilibrium problem£&^£, on both sides of 
the transition, appears to be washed out by nonequilib¬ 
rium effects. Indeed, as we have discussed earlier, the 
scaling of the current in the steady state does not show 
signature of non-trivial power laws, even in the small 
quench amplitude regime where transport is essentially 
set by the quench-induced decoherence scale 7 *, with very 
weak additional dependence from the Luttinger param¬ 
eters. This suggests that the analogy between this scale 
and an effective temperature, meant here as an infrared 
cutoff, cannot be pushed too far. 

Finally it is worth comparing our results for the steady 
state impurit y p hysics with the results and the predic¬ 
tions of Ref. [38l Here the authors computed the zero 
frequency charge susceptibility in a uniform LL after an 
interaction quench starting from a free system. For repul¬ 
sive interactions they found it to diverge at the backscat¬ 
tering wave vector and concluded that, as in equilibrium, 
even a weak single impurity strongly disturbs the homo¬ 
geneous LL or, in an RG language, the perfect chain fixed 
point is unstable. Then, they computed the local density 
of states close to an open boundary and found that it sat¬ 
urates to a non-zero value, with power-law corrections at 
low frequency. From the scaling of these corrections they 
conclude that the steady state analog of the open-chain 
fixed point is stable, with a power-law finite temperature 
conductance, which eventually crosses over to a non-zero 
value at low enough temperature. Based on their results 
the authors argue that the fixed point structure of a single 
impurity in a nonequilibrium steady state LL is similar 
to the one in equilibrium^. 

As we have seen, our transport results combined with 
our analysis of the transient orthogonality catastrophe 
problem^, suggest a rather different picture and high¬ 
light the importance of inelastic effects for the steady 
state impurity physics after the quench. In particular 
in the back-scattering case we found that, due to the 
global quench, a local non-linearity generates an effective 
decoherence in the problem, even in the regime where 
it would nominally be irrelevant at the one loop level. 
This quench-induced decoherence energy scale reduces 
the value of the zero temperature conductance away from 
the perfect uniform limit and, more importantly, com¬ 
petes with the growth of the backscattering as one enters 
the strong coupling phase, corresponding to K < 1 for 
the case Kq = 1 considered in Ref. 0 The result of this 
competition is clearly a subtle issue to establish firmly. 
Our results for the Loschmidt echo^fi suggested that well 
within the putative strong coupling phase this new en¬ 
ergy scale acts as an effective infrared cut off, resulting 
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in an exponential decay of the echo. This result would 
suggest that, quite differently from the equilibrium case, 
a weak single impurity does not disturb substantially the 
uniform system, i.e., one should not trust the perturba¬ 
tive break down for K neq < 1. 

To further confirm this intuition, we have considered 
the opposite limit of an open chain, perturbed by an in¬ 
teraction quench and by a sudden switching on of a local 
tunneling. For any finite bulk quench amplitude we find a 
non-vanishing zero temperature conductance essentially 
set by the quench-induced decoherence scale. This result 
suggests that, differently from the equilibrium case, one 
never reaches the open-chain fixed point when the bulk 
has been quenched because the weak tunneling is always 
effective in generating a non-zero conductance. We may 
interpret this in an RG picture for the tunneling problem, 
dual to what we discussed in Ref. [36| for the backscatter- 
ing case, where the local non-linearity generates an effec¬ 
tive decoherence that cuts off the RG flow. While a non¬ 
zero tunneling conductance was also found in Ref. dU, 
the physics behind this effect was not discussed. We now 
fully clarify its origin as a quench-induced decoherence 
phenomenon. 


V. CONCLUSIONS 

In this paper we have studied the problem of transport 
through a localized impurity in a Luttinger model far 
from equilibrium due to a bulk interaction quench. This 
model in its ground state has a familiar and well studied 
zero temperature transition between a conducting and 
an insulating phase, depending on the value of the bulk 
interaction parameter K , and has an associated finite- 
temperature and finite-voltage correction which shows 
universal power law behavior. 

Here we have presented perturbative calculations in 
the strength of the impurity potential starting both 
from the limit of a uniform liquid and from the one 
of two decoupled semi-infinite systems. In the former 
case a standard bosonization approach gives rise to a lo¬ 
cal backscattering term (described by a boundary sine- 
Gordon problem in a time-dependent nonequilibrium 
bath) that we treat to leading order in perturbation 
theory. In the latter, we use open boundary bosoniza¬ 
tion to treat the problem and obtain a dual formula¬ 
tion again in terms of a boundary sine-Gordon model 
where the non-linearity comes from a local tunneling be¬ 
tween the isolated wires. Our results quite generically 
reveal that the nonequilibrium excitation of bulk modes 
induced by the global quench has important and pecu¬ 
liar effects on the conducting-insulating transition, which 
gets smeared out into a crossover for any small quench 
amplitude K — K 0 7 ^ 0. In addition, the dynamics of the 
current displays a faster decay towards the steady state 
as compared to the equilibrium zero temperature case. 

All together this suggests that the global quench ef¬ 
fectively induces a decoherence mechanism for the local 


degrees of freedom and we have identified an energy scale 
7 * associated with this. While this behavior is qualita¬ 
tively similar to a finite effective temperature, the deco¬ 
herence energy scales enters the problem differently than 
a temperature. In particular the steady state current in 
both limits is set by the decoherence scale but does not 
show power laws in this energy scale. 

An interesting question concerns the generalities of 
these results beyond the quenched Luttinger model. 
We may speculate that quench-induced decoherence is 
a generic phenomenon that occurs in other interacting 
quantum impurity problems coupled to environments 
which are non-thermal at long times but rather flow to a 
Generalized Gibbs ensemble steady state. We leave the 
investigation of this intriguing question to future work. 
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Appendix A: Green’s Functions — Definitions and 
Identities 

We define the contour ordered Green’s function for a 
real bosonic field cf)(t) in the a, /3 = ± Keldysh basis 

G“*(l,2) =-*<^(1)^(2)} (Al) 

where we have, following the convention of Ref. |49| 

G— (1,2) = G T (1,2) = —i{T(j>(l)(j>{2)) (A2) 

G-+(l, 2 ) = G<(1, 2 ) = -i( mm) (A3) 

G+-(l, 2 ) = G>(1, 2 ) = -*< mm) (A4) 

G++(l, 2 ) = G f (l, 2 ) = -i(fmm) (A5) 

We then define the retarded, advanced and Keldysh com¬ 
ponents as 

G fl (l, 2 ) = —i9(l — 2 )([</>(l), </>( 2 )]) (A 6 ) 

G A (1, 2 ) = i0{2 — 1 )([</>( 1 ), </>( 2 )]) (A7) 

G K {m = -i{{m,m}) m 

and find the following useful relation 

2G a p(l, 2) = G K { 1, 2 ) - aG A { 1, 2 ) - /?G fl ( 1, 2 ) (A9) 
from which the following relations follow 

G R ( l, 2 ) = -i^/3G Q/ 3 (l,2) (A 10 ) 

o(3 

G K (l 1 2)-f3G R (l,2) = J2 Ga,»(l,2) (All) 

a 

G k (1 1 2) = IJ2 Ga/j(l,2) (A 12 ) 

a (3 
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Appendix B: Forward Scattering Contribution to 
the Current 

In this appendix we show that the forward scattering 
does not contribute to the time dependent current. To see 
this we follow standard step s 51 ’ 52 and introduce first the 
even and odd combinations of the LL fields 9, 0, defined 
as 


+92 [ dxp L (x)p R (x ) (Cl) 

Jo 

where H {iee = v F J2 a = L ,R s <* fo dxi>t t id x i>a with s L = 
—s R = 1. The crucial observation is that imposing 
open boundary conditions introduces a constraint be¬ 
tween right and left moving fermions which are no longer 
independent, but rather satisfy^ 


0e/o 

^e/o 


(f)(x) ± <f>{—x) 

V2 

6(x) ± 6{—x) 


(Bl) 

(B2) 


Then, it is convenient to introduce new chiral bosonic 
fields ® s / a {x) = y/K9 0 / e (x) + -+=(f) e/o (x) satisfying 
[$ s (a;),$ s (y)] = [$ a (a;), $ a (y)] = ~insign(x ~ y) and 
[4> s (a;), tE» a (?/)] = 0. In terms of these fields the bulk LL 
Hamiltonian decouples in the two channels 


n = ^ I dx Y (^) 2 = Y ^ 


(B3) 


It is particularly convenient to write the local scattering 
in terms of these fields, due to their well defined proper¬ 
ties under inversion. We have 


0 l(~x ) = ~tf> R (x) (C2) 


and is a direct consequence of the existence of a single 
Fermi point. We can therefore fold the line into (—L,L) 
and use Eq. (1C2I) to express the Hamiltonian only in 
terms of right movers 


H 0 = Hfree+y / ^ Pr{ x )Pr( X )+'^ 


dxp R (x)p R {-x ) 


(C3) 

where H( lee = v R f dxip^idxi/jn and one has to notice the 
non-local g® interaction that results from the single chiral 
field representation. The right moving fermions now sat¬ 
isfy periodic boundary conditions and can be bosonized 
as usual in terms of a single bosonic field 


i/j{x) 



(C4) 


Vloc — 9fs \J 2 ^ x < ^a(*^)|a;=0 T ^bs COS a/2A $ s (x — 0) 

i.e. the forward scattering couples to the antisymmetric 
mode while the backward term only to the symmetric 
one. Now we recall Eq. © of the main text, which 
shows that the time-dependent current I(t) is expressed 
only in terms of the Green’s function of the local field 
0(0), which turns out to be related only to the symmet¬ 
ric combination 0(0) = i J A'/2d> s (0), from which we con¬ 
clude that the current does not depend on the forward 
scattering potential. 


in terms of which the Hamiltonian reads 


H 0 = 


Vf 


1 + 


94 


27 tv f 
g 2 _ f dx 
47T / 27T 


dx 

2tt 


(■ 9 x ipi(x )) 2 + 


d x ipi(x)d x (pi(-x) (C5) 


To diagonalize this problem it is useful to decompose the 
field into normal modes 


<p(x) = Y \rr e ~ aq/2 eiqx + b l e ~ iqx ) ( C6 ) 


q>0 


in terms of which the Hamiltonian becomes 


Appendix C: Quenches in the Open-Boundary 
Tomonaga Luttinger Model 

In this appendix we briefly discuss the quench prob¬ 
lem in a TL model with open boundary conditions, and 
compute in particular the local bosonic correlator rele¬ 
vant for obtaining the tunneling current. We first re¬ 
call how to diagonalize the initial Hamiltonian Ho, as 
the same strategy will be used for the Hamiltonian after 
the quench, and this will highlight the major differences 
with the bulk case. We follow the treatment in Ref. [5Cj 
to which we refer the reader to for further details. 

Let us start with the problem of spinless fermions tp(x) 
defined on a line [0,L] with open boundary conditions 
(OBC) with the Hamiltonian 


Ho = Y A °<1 b l h 0 - ( b l b l + b A) (C7) 

q> 0 

where 

A> = v F (l + gl/2Fv F ) (C8) 

= i ( C9 > 

We can diagonalize the Hamiltonian via a Bogoliubov 
rotation parameterized by the angle ao- 

b q = cosh ao g q — sinh ao rj q (CIO) 

b J = — sinh ao r] q + cosh ao r)\ (Cll) 

The Hamiltonian when written in terms of the fields r] q 
is, 


H 0 = H{„ 


+ y / dx (pl(x) + p%(x)) 


Ho = Y v o q vlvq ( C12 ) 

Q 
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for an angle ao given by 


tanh2ao = — 


2B 0 

Aq 


g%/2irv F 
l + g°/2nv F 


where (assuming as usual a Galilean invariance preserv¬ 
ing quench such that vqKq = vK) 

(C13) fq(t) = coshao cose q t — z cosh (2a — ao) sine g t 

g q (t ) =sinhci!o cose q t + i sinh (2a — ao) sin£ g t 


The sound velocity reads 

no = A 0 cosh 2ao + 2 B 0 sinh a 0 = 
= v F \J (1 + gl/2Trv F ) 2 - (g^/2nv F ) 2 

while the angle ao is defined as 


(C14) 


1 + g^/2nv F — go/2nv F ,„ lr \ 

= TW^7TIJ^ sK ° (C15) 


In particular we can write down the time-evolved local 
bosonic field ip(x = 0, t) as 

V>{t) = E \[~T e ~ aq/2 (C23) 

q V y 


where 


B q (t) = 


1 


K 0 . 

cos e a t — sm e a t 


(C24) 


J q\ u ) ~ \fK~ 1 " IF allic -q u 

from which all the local Green’s function can be evalu¬ 
ated. We start with the retarded component 


We can now proceed to discuss the problem of an in¬ 
teraction quench of the bulk interactions g^g® —t ff 2 , 34 - 
This amounts to suddenly change A 0 ,B o —> A,B in the 
bosonic Hamiltonian. We can repeat the previous steps 
and introduce operators 7 ], 7 9 

b q = cosh ao 7 g — sinh ao 7 ^ (C16) 

b\ = - sinh a 0 j q + cosh a 0 7 ^ (Cl7) 

giving a quadratic Hamiltonian when written in terms of 
the fields j q 


G« v (t, t') = - W{t - t'){ W(t)Mt')]) (C25) 

Inserting the time-dependent field expansion (1C23I) . we 
obtain 

Ojr _ p~aq 

O* = ImB,(f)5*(t') (C26) 

q> 0 ^ 

Similarly for the Keldysh component 

G« ip (t,t') = -i(W(t),<p(t’)}) (C27) 

we obtain 


H = E vq 7<? 
q 


for an angle a given by 


tanh 2a = — 


2 B 
~A 


g 2 / 2ttv f 
1 + gi/2iTv F 


such that 


K = exp 2a 


(CIS) 


(C19) 


(C20) 


Using the above transformations, we can easily express 
the time evolved bosonic operators in terms of those di¬ 
agonalizing the initial Hamiltonian 

b q {t) = e lHt b q e~ lHt = f q (t) g q - g q (t) <q\ (C21) 
b\(t) = e* Ht b\ e~ iHt = —g* (t) Vq + f*(t) r,\ (C22) 

I 


G^(t,t’) = J2 — coth (^) ReB q (t)B* q (t’) 

(C28) 

After some algebra and using the fact that 

R eB q {t)B*(t') = — cos e q t cos e q t' + 

K 0 

sm£ q ts'me q t' (C29) 

lmB q (t) B*(t') = — — sin e q (t — t') (C30) 

H A 

we obtain the following results for the local Green’s func¬ 
tions 

2 7r p~ a q 

G*q,{t > t') = —— sin E q {t - t') = 

^ q>0 ^ 

2 

= —— arctan (A (t — t')) 

K 

and 


Gp V (t,t') = — 


-aq F 

K 0 L^ — C ° th( # ) 


E 


COS Eq (t t') ^1 + +COS £q{t + t') ^1 - 


(C31) 
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i 

'^o 


7 e -’^> 


COS £g(t - t') ^1 + + COS £q{t + t') ^1 - 


K 2 

i_£. 

K 2 


(C32) 


The overall result is that the local Green’s functions for 
the boundary field in the OBC case can be obtained from 
the bulk one by sending 


K 0 ,K-> 


2 2 
K~o'K 


(C33) 


where we have defined $ ob (fi,t 2 ) as the following com¬ 
bination of local Green’s functions 


^ 2 ) 


+* \ a G^ v (ti, f 2 ) + bG^ v (ti, f 2 )] 


a result which naturally generalizes the well know equi¬ 
librium result. A simple interpretation^ is that the field 
in the OBC case corresponds to the 9 field in the trans¬ 
lational invariant case (the field cf> being frozen due to 
the boundary condition^), which is known to have cor¬ 
relators with K —y 1/K. On top of this, the field is now 
living in a semi-infinite space, which explains the extra 
factor of 2. Finally, let us compute the following cor¬ 
relator which will be useful in evaluating the tunneling 
current 


fabitiM) = (e -wM) (C34) 

As before, this can be written only in terms of the local 
Green’s function of the field ip a (t) 

f(ti,t.2)ab = exp t 2 )^j (C35) 


Using the above results we can write 

exp [-*T ab (fi,f 2 )] 


fab(t 1 : t 2 ) = 


(l + A 2 (h - t 2 ) 


) IS dual 
neq 


where we have defined 


k: 

as well as 


dual 


neq / tr 4AT n 


K 2 

i±^ 

K 2 


(C36) 


(C37) 


(1 + 4A H\) (1 + 4A 2 t|) 
(1 + A 2 (U + f 2 ) 2 ) 2 


Rtr 1 11 / 2 


(C38) 


T ab (ti,t 2 ) = y {aG^(t 1 ,t 2 ) + bG^(t 1 ,t 2 ))( C39) 
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